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TITLE 

TYPE 

ACCURACY 

DURATION 



DESCRIPTION 



METHOD OF USE 



Automatic Eigenvalue - Eigenvector Program (DOI Only) 
Entire Program 

Depends on the condition of the symmetric matrix. 
Usually 10-11 decimal places. 

(a) Input program M - 3 about 10 seconds 

(b) Input Matrix - about (n + n)/kQ seconds 

(depending on number of digits ) 

3 -3 

(c) Computation time - about ( 5 n ) 10 seconds per 

iteration, -where an iteration means (n - n)/2 
orthogonal transformations . Most matrices require 
from k iterations ( n = 3) to 7 iterations (n = 23) 
for convergence. 
This program computes the eigenvalues and eigenvectors of a 
symmetric matrix of order n < 2k, It is essentially a com- 
bination of library subroutines M-0 (modified), X-l, N-J* R-l 
and P-6. The user of this program need only supply the 
matrix and two parameters . For a description of the math- 
ematical method of computing the eigenvalues and eigen- 
vectors see program M-0. 

(1) Read program tape in until it stops. 

(2) Insert parameter tape in reader and in the two 
parameters . 

(3) Replace program tape in reader and read in until it 
stops a second time. 

(k) Insert matrix tape in reader and read in aaatrix. 

(5) Replace program tape in reader. 
After computation has been completed the remainder of the 
program will be read in and results will be punched out. 
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PREPARATION OF PARAMETER TAPE 

The user of this program must prepare a parameter tape 
as follows: 

00 3K 

00 F 

00 gF 

00 F 
00 nF 

2k 999H 
The parameter g indicates the number of deciaal digits desired in the print 
out and n is the order of the matrix, 
P^PABATION OF MATRIX TAPE 

The symmetric matrix (a. .) must he scaled so that 






ij 



<l/2. 



ij 



Then a. j, for i >J, is punched as a sign digit followed 

by up to eleven decimal digits. Element a is followed by an N* 



EXAMPLE; 



To diagonalize the 3rd order matrix 



= 10" 



1 


<5 


> 


2 


k 


: 


-3 


5 


6, 



.01 


■ i 


,02 


«0+ 


i..03 


,05 



and print answers to 11 decimal places prepare the tapes 



00 


3K I 


00 


F 


00 


11F 


00 F 


00 


3F 


2k 


999^ 



Parameter tape 



Several spaces 

+01 
+02 
+Qk 
-03 
+05 
+06 
n 



Matrix tape 
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Notice the order in which the elements are punched. Since the scale factor 

-2 P 

is 10 the eigenvalues must be "unsealed" by a factor of 10 . Ho scaling 

is necessary for the eigenvectors since they are normalized . 



The Printed Results: 

For the 3^d order matrix above the results are 

+0^00000000^ k iterations 

-O3O773O3608 

+03841390157 eigenvalues 

+10235913^50 

I 

+6527195^698 

-551^5508707 1st eigenvector 

+51947519662 

IT 

+7^65509718^ 

+58486128172 2nd eigenvector 

-31717333989 

N 

-12891407756 

+59^83995162 3rd eigenvector 

+793^3972205 

The first number indicates that k iterations (see above) were required for 
convergence. Ignore the k in the eleventh decimal place produced by round- off , 
The next three numbers are the eigenvalues and following these are the three 
eigenvectors corresponding to the three eigenvalues, in the same order. The 
TX between eigenvectors permits the output tape to be used as an input tape 
for another problem if one desires to do so. 
ARITHMETIC TEST 

In case of an overflow or other type of arithmetic failure the 
program will stop on a 2k 182F order. If started again it will print N(l) 
and N(l4) as 20 consecutive sexadecimal characters and continue * If the first 
ten characters do not agree with the last ten there has been a' failure in 
Arithmetic. 
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TAPE 


MAKE-UP 


1. 


Decimal Order Input No. x-1 




Spaces (parameter tape Inserted here) 


2. 


00 512K 




L5 ^F 




Lk ^F 




ko 5F 




Lk kF. 




LO 12L 




kO 600F ( 




L5 ^F 




Lk 12L 




^0 601F 




50 601F 




i y _rj -' 




S5 F 




10 IF Generate parameters 




Lk 11L 




kO 7F 




50 600F - 




75 k? 




S5 F 




10 IF 




Lk 11L 




^0 13F 




26 999F 




00 F 




00 192F 




00 F 




00 IF 




26 512N 


3* 


00 31K 




Decimal Number Input No, H*2 


K 


00 52K 




Square Root No. R-l 
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5. 00 62K 

Single 'Column Print No. P-6. 

6. 00 20K 
80 40F 
^0 161F 
26 161F 
00 F ' 
00 IF 

00 IF 

00 Sk 

00 Sif 

00 192F 

00 192F Constants 

80 S5 

00 S5 

20 F 

00 F 

OOF 00 0100 

0000 0010 J 

80 F 

00 F 

00 F 

00 F 

JO S7 

Jk ST 

7. oo 76K 

Library Routine M-0 

8. 00 800K 
L5 10F 
LO 1021F 
32 1017F 
LJ 10F 



4»Q* 



9'. 



10, 



22 12QF 

L5 ioeiF 

HO 1QF 
L5 1020F 
HO 2F 

26 122F 
7L 4095? 
LL H095F 
3L kO^F 
LL H095? 
L3 6F 
32 102 3F 

26 111F 

27 113F 
00 850K 
L5 800F 
kO 1015F 
L5 L 

LH 5L 
HO L 
LO Hl 
36 191F 
26 L 
L5 809F 
HO 1023F 
00 IF 
00 IF 
00 &K 
Hi F 
L5 28F 

Ho s? 
L5 9F 
LH 22F 
H6 9F 
LO 19F 
36 lHF 



amendment to No» H2 



Plants above amentlmeni beginning at 2j$l*| 
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50 998P 

24 31F 

00 lUK Generate unit matrix 

L5 P 

L0 18F 

36 8F 

L5 F 

L4 22F 

kO F 

11 63F 

26 9F 

00 S4 

00 S4 

NO F 

L5 9F 

00 191K 

50 18F 
75 18F 
10 IF 
LU 7F 
00 20F 
kS 19F 
>1 191F 
26 8F 

00 998K 

40 192F 
50 5F 
00 6k 

00 F 

2h 76F 

26 85 ON 

Spaces (insert matrix here) 



11. 80028 400K2 
80028 400K3 
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19026 
80028 
L40K2 
80028 
41005 
L5016 
L50SL 
2603F 
LkOO'J 
L0094 
L5007 
40007 
L5094 
L4006 
L1000 
2603F 
L0093 
L5005 
40005 
41005 
L4016 
L0017 
22096 
L5000 
L5000 
92302 
41007 
L4016 
26088 
L0017 
L508K 
92302 
41007 
L501I? 
46094 
00014 
2608O 
2609J 



260K1 
400KL 
400K2 
4007L 

41006 
40007 
50082 

L5082 
46082 
36095 

L4oi6 
26082 
L4005 
46o8k 
5008K 
22098 
3608L 
L4017 
22096 
L5006 
40006 
32092 
OFOOO 
00000 
00000 

92081 
L5007 
40007 
L5007 
3609s 
2608N 
92081 
2608N 
46094 
L500J 
46093 
00000 
00001 



Printing 
Instructions 
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